Black hole entropy and isolated horizons thermodynamics 
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We present a statistical mechanical calculation of the thermodynamical properties of (non rotat- 
ing) isolated horizons. The introduction of Planck scale allows for the definition of an universal 
horizon temperature (independent of the mass of the black hole) and a well-defined notion of energy 
(as measured by suitable local observers) proportional to the horizon area in Planck units. The mi- 
crocanonical and canonical ensembles associated with the system are introduced. Black hole entropy 
and other thermodynamical quantities can be consistently computed in both ensembles and results 
are in agreement with Hawking's semiclassical analysis for all values of the Immirzi parameter. 

PACS numbers: 



Black holes are remarkably simple gravitational sys- 
tems for distant observers so long as one neglects quan- 
tum effects. However, for h ^ their physical behaviour 
remains an open question whose complete answer re- 
quires a full-fledged quantum gravity theory. The most 
difficult challenge is perhaps to unravel the physics close 
to the classical singularity dressed by the event hori- 
zon. The quantum gravity effects are also felt by ob- 
servers outside the event horizon; as clearly indicated by 
Hawking's semiclassical calculations |l| which show that 
a generic BH radiates as a perfect black body at Hawk- 
ing temperature T H proportional to its surface gravity 
and has an entropy S = A/4& where t p = h 1 ' 2 (in units 
G = c = 1) is the Planck length and A is the classical 
area of the event horizon. The analysis of these thermo- 
dynamic aspects of BHs is well within the reach of the 
existing developments in quantum gravity. 

In fact, an account of the thermal properties of BHs 
from the statistical mechanical treatment of the micro- 
scopic germs arising in the underlying quantum theories 
of gravity has now become a standard benchmark for 
testing those theories. In this paper we attack this prob- 
lem from the viewpoint of loop quantum gravity (LQG). 

The problem of computing black hole entropy in the 
framework of LQG has a long history (see [2| and ref- 
erences therein). Despite some small differences among 
various treatments, one common viewpoint has surely 
emerged which is that in order to find agreement with 
Hawking's semiclassical results one must fix the Immirzi 
parameter 7 (a dimcnsionless constant that labels vari- 
ous inequivalent kinematic quantizations of LQG) to a 
critical value 70. Although logically viable, this pecu- 
liar tuning of 7 has arguably become the Achilles' heel 
of the LQG analysis. In this paper we propose an alter- 
native analysis of black hole entropy from LQG whose 
main merit is to reconcile Hawking's semiclassical results 
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with the statistical mechanics treatment of LQG without 
having to fix the Immirzi parameter. 

The key conceptual input is that the first law of black 
hole mechanics needs to be modified from the classical 



form dE^ 



KalA/(8ir) + work terms to 



dE n 



8" 



dA + (AoadN + work terms, 



(1) 



where k is the surface gravity of the event horizon, Eoa 
is the BH-mass measured by the stationary observers at 
infinity 1 and the second term originates naturally from 
the underlying quantum geometry description of the BH 
horizon where the integer N refers to the number of topo- 
logical defects in the quantum isolated horizon (IH) and 
plays the role of a quantum hair for the BH. Then the 
quantity /i plays the role of chemical potential. As can be 
immediately seen, the above modification of the first law 
is fully consistent with standard results for Schwarzschild 
BH if 2 E x = M, ^ = -a( 7 )T 1 co and 



S= \% +N(J{1) 



(2) 



where A is the classical area of IH and (7(7) is some func- 
tion of the Immirzi parameter. In the following we will 
show that the above first law and the entropy @ follow 
directly from the statistical mechanics of the basic quan- 
tum excitations of IH in LQG. The Immirzi parameter 
is completely free and enters the entropy formula only 
through the chemical potential. 



1 In the usual statement of first Law dE^o measures the mass dif- 
ference of two different stationary black holes. Usually, one uses 
the notation Eoo = M where M is the ADM-mass of the BH. 
Here, we use Eoo instead because we will also deal with a local 
version of the first law associated with the stationary observers 
in the interior of the spacetimc for which the energy E will differ 
from Eoo. 

2 In a static spacetime local temperature and chemical potentials 
are obtained from Tlgtt] 1 ' 2 = T"oo and ^gtt] 1 ' 2 = Moo, the local 
temperature is called the Unruh temperature; see Landau and 
Lifshitz, Statistical Physics, Part I, §27. 



In the following discussion, wc use standard coordi- 
nates in which the Schwarzschild metric takes the form 



ds 2 = _ (1 _ 2M )dt 2 + (1 _ ^£)-l dr 2 + ^2 (3) 

r r 

One expects that in the semiclassical limit a spherically 
symmetric quantum black hole of large mass M in a sta- 
tionary state is well-approximated by a Schwarzschild BH 
having a test field in the Hartlc-Hawking vacuum state. 
Equilibrium is sustained by a steady incoming flux of ra- 



diation at Hawking temperature T H 



;/(8ttM) past 



null infinity J"~ and a steady outgoing flux of radiation 
at the same temperature at future null infinity J? + . The 
temperature measured by a stationary observer in the 
interior is the local Unruh temperature 



T(r)=T H (l 



2M, 

r 



-1/2 



(4) 



Classically, this temperature diverges at the BH horizon. 
This is due to the infinite blue-shift of the asymptotic 
energy scales at the horizon. However, in the quantum 
theory there is a universal (independent of mass M) local 
temperature at the horizon. More precisely, for an ob- 
server at r — 2M+e the proper distance from the horizon 
is £ = 2(2Me) 1 / 2 and from ((4]) the local temperature is 



Tr, 



f 2 

2irC 



(5) 



Classically, £ — > as e — > 0. Quantum mechanically, the 
closest proper distance £ must be given by the smallest 
length scale that the quantum geometry can probe and 
hence it must be set by the underlying quantum the- 
ory of gravity. For example, in string theory £ must be 
determined by the string tension a 1 ; in LQG £ ~ j£ p . 
Remarkably, none of the physical result depends on this 
scale; so we do not fix it anywhere. Nevertheless, the 
existence of such a scale makes the local temperature ([5]) 
measured by a stationary observer closest to the BH hori- 
zon universal. This is the relevant temperature for the 
quantum theory of isolated horizon and in its own spirit 
we call it the Unruh temperature. 

The usual global notion of event horizon needs to be 
revised in the context of quantum gravity. The very fact 
that BHs radiate in the semiclassical regime makes the 
definition of event horizon as the boundary of the past 
of future null infinity unphysical; [3[ provides a clear de- 
scription of this viewpoint in a simplified setting. In LQG 
this issue is resolved because one uses isolated horizons 
[J]. IH captures the main physical and local features of 
BH event horizons while being of a quasilocal nature it- 
self. In particular, isolated horizons satisfy a quasilocal 
version of the first law [5| 



dE IH = ^dA- 

Sir 



work terms, 



(6) 



where E IH is a suitable quasilocal energy function and 
k ih is a local notion of IH surface gravity. Neither E IH 



nor k ih arc completely determined in the IH framework. 
More precisely, there are infinitely many possible first 
laws according to the choice of k ih as a function of the 
extensive variables entering the first law which when in- 
tegrated provides a definition of E IH . 

This indeterminacy can only be eliminated by an ap- 
propriate physical input which makes an IH the closest 
representative of a BH spacetimc. In the spherically sym- 
metric case, this input is provided by the Schwarzschild 
spacetime. Indeed, there is a natural quasilocal energy 
that one can associate with the stationary observers in 
Schwarzschild spacetime. The four velocity of such an 
observer is u a = £ a /l£ • S} 1 ^ 2 where £ the timelike killing 
vector field (£ = d t in the coordinate system ((3|)). Then 
the Komar mass integral gives 



E r = — —— I e a bcd^ c u a , 

»7T 
S r 



(7) 



where S r is a spherical section of the r = constant surface. 
The integral gives 



E 






A 



(8) 



when evaluated at r = 1M + e where £ = 1\j2Me. This 
gives a natural notion of energy close to the BH horizon 
(see |6| for a full explanation of why this is the correct 
notion of energy for local observers). Since area is the 
only geometric quantity here, no wonder that the local 
energy is determined by the area where 16n£ provides 
the appropriate scaling. The above analysis provides a 
clear-cut justification for the choice of area in the defi- 
nition of microcanonical ensembles used in [7[ , and more 
profoundly recently in [9(. 

From now on we study the statistical mechanical prop- 
erties of quantum IHs. As follows from the basic LQG 
treatment, we consider a quantum IH to be a gas of its 
topological defects, henceforth called punctures. Using 
(|5J), we take the appropriately scaled IH area spectrum 
as the energy spectrum of the gas. Using the LQG area 
spectrum jlOj 



£|ji,j2---) = (^E\/^ + 1)) 



Jl,j2 



(9) 



where j p taking values from the set {1/2, 1, 3/2, ...} is the 
spin associated with the p-th puncture and we used the 
shorthand notation £ 2 = r )£ 2 

The microcanonical ensemble is defined by an energy 
E = A/(16ir£) where A is the classical area of the IH and 
a number of punctures N. A quantum configuration {sj} 
is given by the number of punctures Sj carrying spin-j 
for all possible values of j. Each configuration must obey 
two constraints 



Ci: ^Vj(iHJ^ = 



A 
8< 



C 2 : J^ Sj = N. 



The number of states rf[{sj}] associated with a configu- 
ration {sj} is 



<*[{*•}] =(£**) ! n 



(2j + lp 






(10) 



We look for the configuration that maximizes the entropy 
log(d[{sj}]) subject to the above two constraints. This 
configuration is obtained from the variational equation 



Slog{d[{ Sj }}) - XSC\ - aSC 2 = 0, 



(11) 



where A, a arc the two Lagrange multipliers. Under Stir- 
ling's approximation, this gives the dominant configura- 
tion 



j- = (2j + l) e - A V^^)- 



(12) 



Summing over all spin values j and using C%, we get 

l = e- <T £(2j + l)e-VjW+ 1 ). (13) 



Denoting by d the value of rf[{sj}] for the dominant con- 
figuration, the entropy S = log d is given by 



S = X 



A 

8^| 



aN where 



a( 7 )-log[^(2 J + l) e - A V7(^ 1 )]. 



(14) 



From (3 = 8S/dE\ N we obtain Lagrange multiplier A as 
a function of (3, namely A = f3£ 2 /(2£). Finally, setting 
T = T v and using (0 we get 

.4 



S= Jp + N( rW> where 



a( 7 ) = log[£(2j + l)e- 27r7 Vj'0-+i)]. 



(15) 



The function 17(7) appear at several places in what fol- 
lows. The chemical potential at T = T v is given by 



T 95 



£ 2 

2ir£ y " 



(16) 



which depends on the fiducial length scale £ and the Im- 
mirzi parameter. 

For further discussion it is instructive to consider the 
same system in the canonical ensemble. The canonical 
partition function is given by 



N\ 



^Ens(^i) 



3s j£, 



(17) 



{sj} 3 



where Ej = £ 2 \J 'j(j + l)/£- A simple calculation gives 
logZ = iVlog[^(2j + l) e -' 3 ^] (18) 



^- log Z at T = T v is a 



and the average energy (E) — 

function of N only; this relates the number of punctures 

to the area 



N = - 



.4 



At* a' (7) ' 



(19) 



Note that for all values of 7 the number of punctures 



0< N < 



Moreover, for a fixed macroscopic area 



A, the number of punctures grows without limit as 7 — > 
while it goes to zero as 7 — > 00. For the entropy we get 

S = -^(-^Z) = lo g Z + pA- f (20) 

At T = T v , this expression is identical to the micro- 
canonical entropy (fT5"j) . Basic formulae allow for the cal- 
culation of the energy fluctuations which at the Unruh 
temperature are such that (AE) 2 /(E) 2 = ff{l/N). The 
specific heat at T v is C — 2N£j£~ 2 a"( / j) which is posi- 
tive. This implies that as a thermodynamic system the 
IH is locally stable. The specific heat tends to zero in the 
large 7 limit for fixed N and diverges as h — > 0. 

We now conclude with some discussion of our results. 

Hawking radiation is a global feature of spacetimes in 
which matter undergoes gravitational collapse and settles 
down to some (scmiclassical) stationary BH geometry. 
Quantum isolated horizons are ignorant about the geom- 
etry outside the horizon and hence are not expected to re- 
produce the thermo dynamical properties of a global BH 
spacetime without additional inputs. Ideally, in LQG cal- 
culations, one should put by-hand the information of the 
scmiclassical quantum states approximating the global 
BH geometry outside the IH as a physical input (such 
states are expected to exist in the large BH mass limit). 
Here we brought in some scmiclassical inputs to the sta- 
tistical treatment of quantum IH by setting the temper- 
ature of the IH at the appropriate blue-shifted Hawking 
temperature ((3J) and using the appropriate quasilocal en- 
ergy ((SJ). When these ingredients are incorporated into 
the entropy calculations, the consequences are striking. 

In addition to the total area of the horizon, the total 
number of punctures (topological defects on the mem- 
brane where surface degrees of freedom live) is also a 
Dirac observable, which must play a role in the statisti- 
cal description of the horizon. 

Our result (20) for the entropy S = /3E + aN is fully 
compatible with semiclassical result of Bekenstein and 
Hawking even for a continuous range of the Immirzi pa- 
rameter. This follows form the fact that (dS/dE)^ = /3 
must be the inverse temperature (|SJ) of the horizon and 
(dS/dN)E = er must be related to the chemical potential 
of the horizon. These suggest that the correct first law 
of a quantum IH mechanics should be dE = TdS + fidN, 
where the new term comes from the quantum hair N 
that arises from the underlying quantum geometry of 
IH, or more precisely from LQG. While comparing our 
entropy with the semiclassical entropy one should note 
that that the latter is inferred from an assumed form 



of the first law, i.e. the entropy S has to be such that 
(35'/9£ , co )... = /^oo holds, where the dots refer to other 
possible macroscopic variables that must be kept fixed. 
Once one equates (3^ = /3 Ha wkin g and E^ = -Eadm, one 
gets the desired expression S = A/A(l. Our entropy fully 
complies with this analysis. The hair TV, whose origin is 
purely quantum geometry, is held fixed in the scmiclassi- 
cal analysis. Hence, the term a(j)N is only an additive 
constant to the entropy and at the semiclassical level, our 
entropy is the same as the one of Bekenstein. This closes 
the gap between the semiclassical analysis and the one of 
the statistical mechanics of IH. 

Our chemical potential /i = —Ta("f) is negative for 
small values of 7 < 70. So long as /i ^ 0, we can lower 
the energy of the horizon at some fixed entropy by adding 
more punctures. That means, large number of punctures 
is favoured. Also, for some fixed energy, the entropy 
maximizes for maximum number of punctures. So large 
number of punctures is also favoured entropycally. This 
shows that TV 3> 1 is the right semiclassical limit of ge- 
ometry. Close to the value 70 of the Immirzi parameter, 
the chemical potential tends to zero and for larger val- 
ues, it becomes positive. For 7 > 70, a quantum theory 
may very well exist mathematically, but it seems not to 
exhibit the right semiclassical behaviour. 

The hair TV has its origin in the underlying quantum 
geometry and hence, the first law of classical isolated 
horizons do not possess this term. Classically, the only 
natural value of the chemical potential is zero, which im- 
plies 1 = 2(2j + 1) cxp(-27T7^j(j + 1)). This fixes the 
value of the Immirzi parameter reported earlier and from 
(20) the entropy S = A/4£ 2 p . This result (with some 
mild differences depending up on the IH model) was ob- 
tained in all previous counting [8[. Our present result 
can clearly reproduce these earlier results. However, it 
differs in many important ways from the existing view- 
point. First of all, in (20) the Immirzi parameter does 
not appear as a multiplicative constant. It appears in an 
additive correction to the semiclassical expression. This 
additive term is the quantum correction to the semiclas- 
sical entropy induced by the quantum hair TV. This result 
is more robust in the sense that the semiclassical results 
are reproduced even when 7 does not exactly obey the 
constraint and the chemical potential is not exactly zero. 
Even to reproduce all earlier results one only requires 
the chemical potential to be only close to zero; more pre- 
cisely TV — > 00 and a ~ ^(1/JV), so that the quantum 
correction to the entropy aN ~ ^(1)- 

The quantum statistical mechanics of isolated hori- 
zon is independent of the ensembles (we have shown 
the equivalence of microcanonical and canonical ensem- 
bles). This is an important characteristic of a statistical 
system in thermal equilibrium when some appropriate 
thermodynamic limit is taken. In general, in absence of 
such a well-defined limit in gravity, one expects that a 
black hole as a statistical system may exhibit features 
that are ensemble dependent. Moreover, the thermo- 
dynamic description is ill-defined because, for example 



for Schwarzschiled black hole, one finds the specific heat 
is negative. For quantum isolated horizons, as we have 
shown here, nothing is needed to overcome these difficul- 
ties. The specific heat is positive and the system is in 
thermal equilibrium. This is the main reason why we be- 
lieve that this is the correct statistical description of IH. 
The limit N — > 00 plays the role of the thermodynamic 
limit in our case (in other words, the semiclassical limit 
and the thermodynamic limit are the same). 

Often the grand canonical ensembles provide more in- 
sights into the problem whose partition function is 3? = 
Y^ z N Z(/3, TV) where Z([3, TV) is the canonical partition 
function and z = exp(/3/i). It is not difficult to see that 
iT = [1 - zfiT)}- 1 where f(T) = £(2j + 1) exp (-/?£,-)■ 
The average energy (E) and the average number of punc- 
tures (TV) are related in the same way as (fT9|) . They also 
show that zf(T) = 1 - ff(l/N), so in the large TV limit 
and for T — T L , the chemical potential is the same as 
(fTS]). The entropy S = (3(E) - (TV) log z + log J° is 



O iJn 



0Q.ogN). 



(21) 



Again, the deviation from the microcanonical entropy is 
small in the large TV limit. However, it is important to 
note that the fluctuations in TV, and hence also in E, 
are ^(1). This signals to the fact that the system is in 
a phase transition region (see for example Pathria, Sta- 
tistical Mechanics, 2nd Ed, §4.5). This phase transition 
must have important significance in the quantum geom- 
etry description of IH. It suggests that a quantum IH ex- 
hibits critical behaviour. This might imply that the semi- 
classical limit of IH is critical (similar to the continuum 
limit in lattice gauge theories where correlation lengths 
diverge). A different theoretical possibility is that this is 
relevant for situations when the IH is placed in the envi- 
ronment of other IHs with which it can exchange topolog- 
ical defects. Such environments arise during black hole 
mergers or quantum mechanical pair productions. The 
behaviour of IH differs significantly from the microcanon- 
ical or canonical descriptions in such situations. We keep 
this important issue of phase transition for future study. 

Since the gas of punctures is in equilibrium at a high 
temperature (of the order of Planck mass), statistically 
the punctures may very well be bosonic or anyonic (the 
departure from Boltzmann statistics should be small). 
This may have important implications at other temper- 
atures, especially when the semiclassical approximation 
breaks down and we have to deal with the quantum me- 
chanics of punctures directly. We wish to analyse this 
aspect of the gas also in the future. 

The quantum corrections found here is different from 
the log-corrections that arise in all counting (the latter 
corrections arise from counting the sub dominant config- 
urations and are also present in our model) . A somewhat 
similar log-correction arise in the case of grand canonical 
ensemble and in that case the two corrections compete 
with each other. Investigations of this aspect is again 
kept for a future study. 
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